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I. INTRODUCTION 

The low-energy dynamics of the process 77r — > tttt are 
governed by the Wess-Zumino-Witten anomaly At 
leading order in the chiral expansion the amplitude is 
fully determined by the electric charge e, the pion decay 
constant = (92.21 ± 0.14) MeV [3], and the number of 
colors TVc 



eN 

= (9-78 ± 0.05) GeV'^ 



(1) 



While the analogous low-energy theorem for the neutral 
pion decay 7r° — )• 77 Q, 



" 12^' 



(2) 



concords with experiment to a remarkable accuracy (see 
Ref. [1] for a recent review), agreement between the 
chiral prediction ([T|) and experimental data has long 
proven elusive. Information on i^3^ can be extracted 
from a Primakoff reaction where a charged pion scat- 
ters off the Coulomb field of a heavy nucleus, giving 
access to the cross section for jtt^ 7r~7r" and thus 
to the chiral anomaly. Assuming that the amplitude 
were solely given by Eq. ([T]), the Primakoff experiment 
in the threshold region performed at Serpukhov [6| led to 
— (12.9 ± 0.9 ± 0.5) GeV~^, suggesting some tension 
with the low-energy theorem. 

Subsequently, it was shown that neither the one- 
loop and two-loop [s^ corrections, nor the inclusion 
of vector mesons can fully resolve this discrepancy, 
while a treatment based on dispersion relations was found 
to suffer from the occurrence of unknown free param- 
eters [l3|- A study of leading- logarithm contributions 
shows negligible corrections beyond one loop [ll|. For 



•Electronic address: 'hoferichter@itp. unibe.ch] 
tElectronic address: kubis@hiskp.un i-bonn.de| 
■l-Electronic address: s akkas®hiskp.uni-bonn.de| 



theoretical approaches based on quark-loop and Dyson- 
Schwinger-type calculations, see Ref. and references 
therein. Finally, large electromagnetic corrections in 
77r~ — > 7r~7r° predominantly due to the t-channel ex- 
change of a virtual photon were discovered in Ref. (13| . 
which in combination with higher-order chiral corrections 
reduce the value of the inferred chiral anomaly to F^tt — 
(10.7 ± 1.2) GeV~^ and thus remove the tension with the 
low-energy theorem. Later on, the prediction for the chi- 
ral anomaly was also tested in vr^e^ — > 7r^e~7r'^ [14, 
leading to an extracted value of F^-^ = (9.6± 1.1) GeV~ . 
Both results are now compatible with the theoretical pre- 
diction, which, however, is only tested at the 10% level. 
Therefore, it would clearly be desirable to improve this 
accuracy and confront the low-energy theorem using bet- 
ter data and a refined theoretical approach. 

Presently, a Primakoff measurement of jTr~ — > tt^tt" is 
being analyzed at COMPASS, with the aim of determin- 
ing the cross section not only in the low-energy regime 
as in the Serpukhov experiment 0, but also for center- 
of-mass energies well beyond, in particular including the 
peak generated by the p(770) resonance However, 
the chiral amplitude will only be valid in the low-energy 
region, so that most of the data would be lost for the 
anomaly extraction. In this paper we propose a disper- 
sive framework that can be used to deduce the chiral 
anomaly from a fit to the full data set up to roughly 
1 GeV, incorporating the physics of the p{770) by means 
of the TTTT P-wave phase shift, and thus vastly improv- 
ing the statistical accuracy of the resulting anomaly de- 
termination. Moreover, we carefully revisit the assump- 
tions required to perform the extrapolation to the chiral 
limit and study the internal consistency of our disper- 
sive framework by comparing representations with a fi- 
nite and an infinite matching point. 

Apart from the test of chiral dynamics renewed interest 
in 77r — > tttt has been triggered recently by its impact on 
hadronic light-by-light scattering, whose uncertainty may 
soon dominate the error budget in the theoretical predic- 
tion of the anomalous magnetic moment of the muon (see 
Ref. [l^ for a review). More precisely, one of the most 
important contributions, the ir^ pole term, will be de- 
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termined by the doubly-virtual decay tt" — )■ 7*7* and 
thus the corresponding form factor F^o^.^. (M^o, Qi, (zl)' 
where q^^^ refer to the photon virtualities. For one of 
the photons being on-shell, i.e., the process 7r° 77*, 
one may again write down a dispersion relation based on 
the two-pion cut, and the 777 — >■ tttt amplitude serves as 
vital input for such a representation. (See Ref. fl?! for 
the same arguments in the context of the w, — > i:^^* 
transition form factors.) 

This article is organized as follows: in Sec. HH we first 
introduce the notation and conventions necessary for the 
dispersive representations for 77r tttt that we derive 
and compare in Sec. lIIII The consequences for the extrac- 
tion of the chiral anomaly and the relation to the muon 
anomalous magnetic moment are discussed in Sees. IIVI 
andfVl before we close with a summary in Sec. IVIl 



II. 



KINEMATICS AND PARTIAL- WAVE 
DECOMPOSITION 



We decompose the amplitude for the process 

7(9)71"" (pi) tt^(p2)tt°{po) (3) 
in terms of the scalar function J-{s, t, u) according to 

M{s,t,u) = ie^,^ap^'"PiP2Po^{s, t, u) (4) 

and Mandelstam variables chosen as s = {q + Pi)^, t = 
(pi— ^2)^7 and u = (pi—pq)'^. Working in the isospin limit 
with = , we have on the mass shell s + t + u = 
3Af^, and in the center-of-mass frame we may write 



t = as + bsZ, as - bgZ, 



(5) 



3A1 



,f 2 



Oc — (T„ , 



a: = \ i- 



with scattering angle z = cos 6. By virtue of isospin sym- 
metry, the scalar function J-{s,t,u) is fully symmetric 
in its arguments. The partial-wave decomposition then 
takes the form [ISj 



J-(s,i,u)= fi{s)Pl{zl 

odd I 



(6) 



where Pl{z) denotes the derivative of the Legendre poly- 
nomials. In our analysis we will only be concerned with 
the P-wave /i(s), whose projection formula reads 



his) 



dz(l - z^)F{s,t,u). 



(7) 



In the absence of inelastic contributions, its imaginary 
part is given by 

Im/i(s) = a:{t\{s)yh{s)e{s - 4Af^), (8) 

with the TTTT P-wave amplitude 



t\{s) 



(9) 



parameterized in terms of the phase shift Sl{s). The uni- 
tarity relation (IST) immediately implies Watson's final- 
state theorem [1^], namely that the phase of /i(s) coin- 
cides with ^i(s). Finally, the formula for the cross section 

"^^^ = Ma^s J )l-^(-.M)l 

(10) 

in the isospin limit should be augmented by 



2e2p2 
T{s,t,u) ^Ps^ 



t 



(11) 



in order _to include the dominant electromagnetic cor- 
rection [13]. Strictly speaking, this correction depends 
on the anomalous tt" — 77* form factor as well as the 
pion vector form factor F^{t). However, higher-order 
terms in the radiative corrections will be of minor impor- 
tance, so that the leading form of the chiral prediction, 
see Eqs. ([1]) and (0), which has been employed in Eq. (|lip . 
should prove adequate. Likewise, subleading electromag- 
netic corrections were shown to be numerically irrelevant 
in Ref. ^} 



III. DISPERSIVE REPRESENTATIONS 

In any application of dispersion relations one of the 
basic assumptions concerns the behavior of the spectral 
function at high energies, i.e., the convergence properties 
of the dispersive integral. As elastic unitarity only de- 
termines the imaginary part below the onset of inelastic 
channels, it is crucial that the uncertainties associated 
with the high-energy input be carefully investigated. In 
the infinite-matching-point setup these uncertainties are 
reflected in the behavior of the phase shifts at high en- 
ergies as well as the contributions from inelastic chan- 
nels. It then needs to be shown that the final results are 
insensitive to these assumptions, which will be fulfilled 
provided the integral converges sufficiently fast. In con- 
trast, in the finite-matching-point formulation the phase 
shifts are needed only in a finite energy domain, while 
the imaginary part above the so-called matching point 
Sin is taken as input. The uncertainties generated by 
the input for the spectral function above Sm correspond 
to the uncertainties induced by the high-energy region 
in the dispersive integral in the infinite-matching-point 
case. In this section we will derive and compare both 
approaches for 771 — >■ tttt. 



^ One might be tempted to think that t-channel u) exchange should 
also yield a significant isospin-breaking contribution, given the 
enhanced ujir^'y coupling. In a hypothetical measurement of 
777" — > n'^iT~ , this would be a very sizable effect, enhanced com- 
pared to the p-uj mixing signal in the pion vector form factor by 
about a factor of 9. However, as a t-channel contribution, it is 
suppressed at the sub-percent level even in the threshold region 
for 77r~ — > iT~n^, and becomes entirely insignificant as soon as 
the energy increases. 
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A. Infinite matciiing point and angular averages 

Neglecting the imaginary parts of partial waves with 
angular momentum / > 3 the amplitude J^(s, t, u) may 
be decomposed as (cf. Ref. |20l]) 



T{s, t, u) = F{s) + F{t) + F{u). 



(12) 



In the spirit of the "reconstruction theorem" in the con- 
text of chiral perturbation theory (ChPT) [21], it can 
be shown that this decomposition holds exactly up to 
corrections of chiral power 0{p^'^). The generalization 
of this representation including the absorptive part of 
the F-wave can be found in Ref. [2^, but in this article 
Eq. ([T^ will be sufficient. F{s, t, u) fulfills the once- and 
twice-subtracted dispersion relations Q 



J^(s, t, u) 



1 



ds' 



t 



TT 74JV/2 S' [_ S' — S S' — t S' — U 



(13) 
Im/i(s'), 



Im/i(s'). 



By virtue of Eq. ([T2|). these dispersion relations corre- 
spond to 



1 



= C2 + - -T 



ds' 



TT 74M2 S'^ !_ S' — S S' — t S' — U 



, , Ci 1 f "^ ds' s 

3 TT J4jvf2 S' S' ~ S 



(14) 



i(cW-KCf.) + ir ^^Im^(.') 



with 



(15) 



(16) 



In addition, we define 

with the hat function 

J-(s) - 3((1 - z^) {z^T} - i dzz"T{t), (17) 

being real on the right-hand cut, so that due to elastic 
unitarity ([8]), 

Im/i(s) = Im J"(s) (18) 

= {Tis)+f{s))0{s - AMI) sinJi(s)e-**'i'(^). 

The solution of this equation for the once- and twice- 
subtracted versions becomes [1^ [1^ 

3 TT Ji^p s'{s' - s)\n{s')\ 



a 



(1) 



(2) 



= f](s)^(l-f7(0),s)+ g 



(19) 



2 roc 



Jam- s'^{s' - s)\n{s') 



with the Omnes function 



r2(s) = exp I — 

[ I" J 4 



ds 



4Mi S'{s' - S) 



(20) 



and its derivative fl(s). In combination with Eq. (|17l) . 
these relations allow for an iterative calculation of J-{s), 
which, in turn, determines /i(s) by means of Eq. (jl6p . 
One may check explicitly that this representation fulfills 
Watson's theorem, and eventually leads to 

\Ms)\^ j^{s) cos Slis) 

+ l"(^)l|T + dM/^ s'is'-s)ms')\\ 

= :f{s) cos 5lis) + \n{s)\l^{i - nio)s) 



c: 



(2) 



r . , F[s')sin5\{s') 



^/4m/" s'^{s' - s)\n(s')\ 



(21) 



where the dash denotes the principal value of the inte- 
gral. The full amplitude can then be reconstructed via 
Eqs. (IT3t . (|T5|) . and Watson's theorem. 



B. Finite matciiing point and kernel functions 

The key point of the framework presented in the pre- 
vious subsection concerns the fact that the calculation of 
the angular averages in Eq. (jl7p has to be done numer- 
ically. In an alternative approach, which proves conve- 
nient for the finite-matching-point formulation, first the 
original dispersive representation (|13p is inserted again 
into Eq. (|17p . the angular integral is performed analyt- 
ically, and only in the last step is the final solution ex- 
pressed in terms of Omnes functions (for details of this 
procedure we refer to Ref. ^24*]), so that in the end the 
angular integration becomes traded off in favor of an- 
other dispersive integral. In this way, the result for the 
modulus, 

|/i(s)| = Ai(s)cos5i(s) 

"''"^^^r^-'dM/^ s'[s'^sms')\ 



TT J S'{S' - S)\n{s')\ 



= A2(s) cosdlis) + \n{s)\\ C2(l - 5i(0)s) 



'If'" d.' ^2(s')sin<5i(.') 



^Un s'^is' ~ s)\nis')\ 



+ 



ds' 



Im/i(s') 



s'2(s'-s)|17(s')| 



(22) 
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indeed resembles Eq. d?!]) . but now the inhomogeneities 
A„(s), n € {1,2}, are not given by an angular integral, 
but by 



1 r°° 

A„(s) = - / ds'i^„(s,s')Im/i(s'), 



(23) 



with kernel functions 
3 



i^2(s,s') = ^l(s,s') + 



3A1 



/f2 



„/2 



(24) 



where 



?o(^) 





da; 




z — a; 




1 + z 


1~ z 



(25) 



denotes the lowest Legendre function of the second kind. 
In addition, we have introduced a finite matching point 
Sm; i.e., Im/i(s) is assumed to be known for s > Sm- The 
Omnes function 



n{s) 




4M2^* s'{s' - s) 



(26) 



now only involves the phase shift below Sm, but exhibits 
a cusp in the vicinity of Sm, 



n{s) 



(27) 



which complicates the numerical implementation, cf. 
Ref. [13]. Once this is accounted for, Eq. ([22]) pro- 
vides an iterative scheme that determines |/i(s)| for 
4M^ < s < Sm, and thus, by means of Eq. (fT3|). the 
full amplitude J-{s,t,u).'^ 

Another interesting feature of the finite-matching- 
point derivation concerns the fact that statements about 
the strict validity of the final integral equation in the 
context of the convergence of the partial-wave expan- 
sion are possible. Given certain analyticity properties 
of J-{s^t,u), one may calculate the corresponding maxi- 
mally allowed value of the matching point; e.g., starting 
from fixed-i dispersion relations and assuming Mandel- 
stam analyticity [25j . we obtain 

W = (1.2GeV)^ (28) 

which already indicates that a dispersive representation 
cannot be rigorously valid far above 1 GeV. Further 
details of the derivation of s^ax are collected in Ap- 
pendix |X1 




FIG. 1: Cross section for -yn rnv for C2 = 9.78 GeV"^ (solid 
line) and C2 = 12.9 GeV"^ (dashed line). The curves for the 
finite/infinite- matching-point setup lie on top of each other. 



C. Comparison 

The free parameters of the dispersive representations 
presented in the previous subsections will be determined 
by fitting to the experimental result for the cross section 
ct(s) once available. With two subtractions, these are the 



constants c!^^^ and 6*2^'' in the infinite-matching-point 
case and C2 for a finite matching point. Moreover, absent 
independent information on the imaginary part, Im/i(s) 
will simply be set equal to zero above Sm,'^ so that Sm can 
be regarded as another free parameter that determines 
at what energy the amplitude is guided to zero. From 
this point of view, both representations involve two free 
parameters that need to be fitted to experiment. 

As a starting point for the numerical implementation 
we take T{s,t,u) ~ 0, but the outcome of the itera- 
tion is completely insensitive to this choice. (We have 
checked that e.g. using a vector-meson-dominance ansatz 
as starting point leads to exactly the same results.) The 
solutions for |/i(s)| and ^(s), for the finite- and infinite- 
matching-point formulations, respectively, are obtained 
by iterating the integral equation until these functions 
change by less than 10~^, which, for the infinite match- 
ing point, requires six iterations. We find that the it- 
erative scheme for the infinite matching point converges 
faster than its finite- matching-point equivalent: the num- 
ber of iterations needed for a given accuracy is reduced 



.(2) 



^ The dispersive treatment in Ref. llOU employs this kernel-based 
approach in the infinite-matching-point case. Since the behavior 
of the inhomogeneities A„(s) for large s induces a rather slow 
convergence of the (second) dispersive integral for s — > 00, we 
only use this formulation in combination with a finite matching 
point. 



^ Indeed, the (preliminary, uncorrected) count rates presented in 
Ref. [Tsl indicate that at high energies the cross section becomes 
essentially zero. If the final cross section at high energies turned 
out to be non-negligible, the corresponding effect could be in- 
cluded by modeling the high-energy tail of /i in such a way that 
the experimental cross section be reproduced. Either way, for a 
sufficiently high matching point the influence on the energy re- 
gion < 1 GeV is expected to be of minor importance, cf. Ref. |24| . 



5 



by roughly a factor 2, and the comphcations due to the 
cusp of the finite-matching-point Omnes function (j27p 
are absent. On the other hand, for the implementa- 
tion of the infinite-matching-point scheme, one needs to 
specify an integration cutoff A and study the sensitivity 
to a particular choice of A. In the following, we take 
A = 1.8 GeV, which proves sufficiently large for the en- 
ergy region < 1 GeV to remain unaffected under cutoff 
variations. 

For illustrative purposes, we now take Sm = (1.2 GeV)^ 

(2) 

and fix in such a way that for a given C2 the infinite- 
matching-point amplitude vanishes at Sm. We include 
the electromagnetic corrections according to Eq. ([TT|). 
although these are only relevant in the low-energy re- 
gion where the cross section is rather small (see Fig. 2 
in Ref. [13); so that these corrections mainly affect the 
conversion between a(s) and the total cross section of 
the Primakoff reaction [13]. Taking Sl{s) from the anal- 
ysis of Roy and Roy-like equations in Refs. [2^ [27| . 
we obtain the cross section depicted in Fig. [1] where, 
as expected given the choice of C2 , the finite/infinite- 
matching-point results are indistinguishable on the scale 
chosen. The only difference concerns the behavior above 
Sni, where the finite-matching-point curve stays at zero, 
whereas the infinite-matching-point cross section slowly 
starts to rise again. This example nicely demonstrates 
the inherent consistency between the two methods, and 
we expect that as soon as data are available, the compar- 
ison between both fits should allow for a valuable consis- 
tency check regarding the extracted value of the chiral 
anomaly. Moreover, Fig. [T] shows that already moderate 
changes in C2 entail large differences in the height of the 
p peak, so we expect an enhanced sensitivity to the chiral 
anomaly if indeed the whole spectrum is included in the 
analysis. 



with 

(30) 

■^ttAm^ s'2(s'-s)|f7(s')| /' 

where ^i(s), T2^\s), and ■p2'^\s) have been rescaled ac- 
cordingly and are all determined from J-i{s), J'2^\s), 
and Jg^' (s) according to the angular averaging prescrip- 
tion (|17p . Obviously, the basis functions as given in 
Eq. pO|) can be calculated iteratively, independently of 
any subtraction constants, so they can easily be pro- 
vided to experimental analyses in numerical form, allow- 
ing for a straightforward extraction of the subtraction 
constants given as linear prefactors in Eq. (|29)) . Fur- 
thermore, Eq. (|29p explains the strong dependence of the 
cross sections shown in Fig. [T] on the subtraction con- 
stants: they really serve as multiplicative constants in the 
whole (elastic) energy range, yielding the desired strong 
impact of the p resonance region on the extraction of the 
anomaly. 

We finally remark that for the finite-matching-point 
scenario Eq. (j22p . a similar linearity property can be 
used as long as the input above the matching point is 
set to zero (as we have done above). For varying in- 
put above the matching point, as well as for different 
matching points, different basis functions would have to 
be calculated. 



IV. EXTRACTING THE CHIRAL ANOMALY 
FROM CROSS-SECTION DATA 



B. Matching to chiral perturbation theory 



The extraction of the chiral anomaly from experimen- 
tal cross-section data will proceed in two steps: first, the 
free parameters of the dispersive representation(s) of the 
previous section, the subtraction constants, need to be 
fitted to data; second, these subtraction constants have 
to be related to the chiral anomaly by matching to ChPT. 
We describe these two steps in the following subsections. 



A. Fitting dispersive representations to data 

For the practical purpose of a fit to data, the rep- 
resentations Uni) can be significantly simplified, based 
on the observation that they are completely linear in 



the subtraction constants C'l and C 
tively We can rewrite 



(1) 
2 ' 



c: 



(2) 



respec- 



(29) 



The chiral anomaly is defined as the value of J-{s, t, u) 
a,t s = t = u — 0. Since the dispersion relations 
have been derived using on-shell kinematics, they cannot 
be used to perform the required extrapolation to the chi- 
ral limit. (For a more detailed discussion of this point 
see Appendix [Cl) Therefore, the appropriate procedure 
to extract the chiral anomaly involves on-shell matching 
to ChPT, with the chiral expansion thereafter used to 
extrapolate to the chiral limit. 

The one-loop, 0(V'), ChPT amplitude may be ex- 
pressed in the form [3| 

J-(<^) (s, i, u) = (s) + J-(s) (<) + ^(6) (u). 



4M2 



,w}V) 

S'2(s'-S) 



(31) 



F:- 



37r< 



647r^ 
~3e" 



967r2F2 



1 -l-loa 



M: 
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where C2(^) is an 0{p^) low-energy constant (LEG), and 
Im f[^\s) denotes the one-loop imaginary part 

Im/f = aU['\s)ti^Hs) = -:^3.9^. (32) 



The dispersive representations ()13p may be summarized 
in the form 



(33) 



^ I f°° ds' ( e . ^ , . 

c + —\z^. + zr-. + zr-M^h{s'). 



with C ~ C2 for two subtractions and 



c = c. + HSf ,34) 

1" JiMl S'^ 



for the once-subtracted case. After the fit of the dis- 
persive amplitude to cross-section data, either in the 
finite- or infinite-matching-point scenario, the partial 
wave /i(s) as well as the constant C will be known. De- 
manding that 

J-°^(s,t,w) = J-(s)(s,i,u) + 0(p*) (35) 

in the low-energy region shows that 

C = F^^{l + ?,MlC), 

since the difference in the integrals is of higher chiral 
order. In consequence, one cannot disentangle the chi- 
ral anomaly from its quark-mass renormalization without 
further input for C2(/i), e.g., its resonance-saturation es- 
timate ^] 



C^2(m) = 



3e 



1287r2M2 ' 



(37) 



where Mp denotes the mass of the p{770), and the renor- 
malization scale is identified a.s fi — Mp* 



Note that the estimate Eq. I I37I I implicitly makes use of the KSFR 
relation |29|| . Allowing for deviations in the relation between the 
pTTTT and the p-y couplings (see, e.g., Ref. [30|) reduces C2{l^) by 
up to 20%, still within an estimated uncertainty for the anomaly 
extraction at the percent level. We are grateful to Stefan Leupold 
for pointing this out to us. Alternative approaches to estimate 
this low-energy constant, based on a chiral quark model [sil or 
a Dyson-Schwinger equation [33 | , seem to us rather inconclusive 
at present. 



Similarly, one finds at two loops Q 

(S, U) = {s) + ^(8) {t) + ^(8) (u) , 



3 i-oo 



ds 



,Im/f)(.s') 



JiMl s'^(s' - s) ' 
^= (4^{60li|+(4^^^^' (^^^ 



where C2 and (J2 are the new LECs at 0{p^). To derive 
the matching condition, we rewrite the dispersive repre- 
sentation as 



J-°^(s, t,u) = C + Csr(s2 + f + u^) 



I r ds 



TT 74M2 S''^ [ S' — S s' — t s' — U 



4AQ ^ 

with the result that 

J-°^(s, t, li) = (s, m) + 

implies 

Fg^l + 3M^C') = C, Fs^D = CsR. 



Ini/i(s'), 
(39) 

(40) 
(41) 



Due to the additional LECs C2 and d2 matching at two- 
loop level does not allow for an unambiguous extraction 
of the chiral anomaly either. However, estimating these 
0{p^) LECs by order-of-magnitude arguments, these ad- 
ditional constraints may prove valuable in assessing the 
systematic uncertainty in the anomaly extraction, in par- 
ticular by comparing the results of the one- and two-loop 
matching procedures. 

To summarize, the chiral anomaly would be derived 
from the dispersive fit using Eqs. ([55]) and ([57]) . po- 
tentially supplemented by the two-loop matching con- 
ditions (|4T|) . According to Eq. ([37)) . C2{fJ.) amounts to a 
4.9% effect, while the chiral logarithm a.t ^ = Mp gener- 
ates another 1.8% correction. In view of the experience in 
the normal-parity sector that resonance saturation works 
rather well for vector mesons |33| . one would expect an 
uncertainty in i^3^ due to C2(/i) at the percent level. 
We stress that this particular source of uncertainty is by 
no means a limitation inherent solely to our approach, 
since any dispersive framework, and of course also a fit 
of the ChPT expression directly to data in the low-energy 
region, will face the difficulty of determining the LECs 
pertinent to the quark-mass renormalization. 
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V. RELATION TO THE MUON ANOMALOUS 
MAGNETIC MOMENT 

The pion transition form factor .F^o^*^. (M^, , gf), 
which determines the strength of the pion-pole contribu- 
tion to hght-by-hght scattering, may be split into com- 
ponents with definite isospin according to [34| 



[Mi ,qi,qi)^ {qi ,qi) + (gi ^ q^), (42) 



where the first/second index refers to isovector [v) and 
isoscalar (s) quantum numbers of the photon with mo- 
mentum qxjqi- If we assume the isovector spectral func- 
tion to be saturated by two-pion intermediate states 
(neglecting heavier contributions from 47r, KK etc.), 
the form factor f^o^{s) — Fys{s,0) fulfills the once- 
subtracted dispersion relation (cf. Ref. (201] ) 



127r2 



ds 



,<z3(,')(F/(s'))V: 



s'3/2(s' - s) 



(43) 

where g^(s) = - and f-j,a^{Q) = Ft,^~J2. In 

this way, the determination of /i (s) from 77r — > tttt cross- 
section data would help in constraining the 7r'^7*7* tran- 
sition form factor for on-shell isoscalar photons. Assum- 
ing an unsubtracted version of Eq. ([^S]) would entail a 
sum rule for /^o^(0), which we discuss in more detail in 
Appendix [B] 

Moreover, F^s fulfills the dispersion relation 

T? ( \ ^ ^ I *2 /"^ , ImF„s(si,x) 
-F-»s(si,S2) = /7r07(si) + — / da; — , (44) 



which would fix the S2 dependence of F^g if the spectral 
function were known, and thus, by means of Eq. (|42l) . the 
whole functional form of Ft^o^*^* . Here, the — >■ tttt in- 
put provides the subtraction function necessary to write 
down a subtracted version of the dispersion relation, and 
therefore it should appreciably improve the convergence 
of the dispersive integral. 



VI. SUMMARY 

In this article, we have shown that a rigorous descrip- 
tion of the anomalous process -> tttt can be extended 
well beyond the chiral regime up to about 1 GeV, using 
methods from dispersion theory. Two different variants 
have been employed and shown to yield equivalent re- 
sults: an infinite-matching-point formulation with two 
subtraction constants, and the use of a finite match- 
ing point, where one subtraction seems to be sufficient. 
Matching to chiral perturbation theory allows us to relate 
the subtraction constants to the chiral anomaly, which 
can therefore be extracted from a fit of the dispersive 
representation to all data below 1 GeV, including in par- 
ticular the large contribution of the p resonance. This 
should allow for a vast statistical improvement of the de- 
termination of the anomaly in the upcoming analysis of 



(I) 



(11) 



FIG. 2: Box graphs constraining the boundaries of the double- 
spectral regions. 



data taken by the COMPASS experiment, using the Pri- 
makoff effect. 

We have pointed out that a reliable determination of 
the 77r — >■ tttt amplitude serves as a vital input in a fu- 
ture dispersive analysis of the 7r° transition form factor, 
which in turn is an important ingredient in the hadronic 
light-by-light-scattering contribution to the anomalous 
magnetic moment of the muon. Given its far-reaching 
impact, we strongly encourage a thorough experimental 
analysis of this process along the lines presented in this 
article. 
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Appendix A: Domain of validity 

In this appendix we briefiy sketch how the convergence 
of partial-wave expansion and projection imposes limits 
on the choice of the matching point. For a more thorough 
derivation we refer to the literature on Roy [s^ [s^ and 
Roy-Steiner 2J| equations. 

Mandelstam analyticity states that the scattering am- 
plitude T(s, t) can be written in terms of double-spectral 
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functions psu, Ptu, and pst as 



1 



+ ■ 



1 



ds'dw' 
dt'dw' 
ds'dt' 



Psu{s',u') 

(s' — s){u' — u) 

Ptujt'jU') 

[f ~t){u' ~u) 

Pst{s',t') 



(Al) 



The integration boundaries follow from the box diagrams 
shown in Fig.[2l which represent the lowest- lying interme- 
diate states possible to contribute to the double-spectral 
functions. More precisely, diagrams (I) and (II) define 
the boundary of the support of pst as the intersection of 

bi{s,t) = 8Af2(2s -I- - t{s - 4Af2) = Q, 
buis,t) = 4M2(s + 2Af2) _ t{s - l&Ml) = 0, (A2) 
or, equivalently, by 

t = Tst (s) = min {Ti (s) , Tii (s) } , (A3) 



where Tj, Tn follow from solving Eq. \K2\ for t. By virtue 
of crossing symmetry, the equations for psu and ptu are 
obtained by the pertinent permutation of Mandelstam 
variables. 

Starting from fixed-t dispersion relations,^ these 
double-spectral regions limit the validity of the result- 
ing integral equation as follows: first, the partial- wave 
expansion of the imaginary part inside the dispersive in- 
tegral converges only for scattering angles z within the 
large Lehmann ellipse [S^l, which is constructed as the 
largest ellipse in the complex z-plane which does not 
reach into the double-spectral regions. In view of Eq. ([S]), 
this automatically constrains the range of allowed values 
for t. Second, a specific value for t is only allowed if the 
corresponding line in the Mandelstam plane avoids the 
double-spectral regions as well. 

The analysis of the second condition leaves the range 



- (17 + l2V2)Ml « -34Mj < t < AM^ 



(A4) 



of allowed values for t (see Fig. [3]), while the Lehmann- 
ellipse constraint can be translated into 



TU{s') = 3M2 - s' - Tst{s') < t < Tst{s') 



(A5) 



for all s' e [4Af^,oo). In fact, this condition amounts 



to exactly the same range already given in Eq. (jA4p (see 
Fig. [3]) . Taking into account that by virtue of Bose sym- 
metry we only need half the angular range to perform 
the partial- wave projection, i.e. t £ [og, Os -I- fos]? we find 
from comparison with Eq. (jA4|) 



smax = (37 + 24V2)m2 « llMl « (1.2 GeV)2. (A6) 



^ It may be possible to further extend the range of vaUdity by con- 
sidering dispersion relations in the manifestly crossing-symmetric 
variables x = st + tu + us and y = stu instead [37{ . 



40 



20 



"5 



-20 



-40 









^\ Psu 












20 30 

s' I Mi 

FIG. 3: Double-spectral regions (top) and allowed range of t 
(bottom) for 77r — >■ tttt. The red lines refer to t = AM^ and 
t = — 34M^, respectively. 



Appendix B: Sum rule for f^o^{0) 

The unsubtracted version of Eq. (j43p implies the sum 
rule 



Ao^(O) 



77 



1 



127r2 



AMI 



(Bl) 

The existence of such a sum rule can be made plausible by 
the following arguments. In the hidden-local-symmetry 
(HLS) formalism [s^, [i^l with the simplest choice for the 
anomalous HLS couplings 03=04 = ! [iO, |4l| the — >■ 
TTTT amplitude reads 

J^{s,t,u) = Fs.jl + ^{Dp{s) + Dp{t) + Dp{u) - 3)|, 
Dp{s) 



(B2) 



Introducing a finite width Fp and approximating the vec- 
tor form factor in a similar fashion, the dominant (pole) 
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part of the integrand in Eq. (jBl[) becomes 



-F, — 

2 (s' 



Mir 



(B3) 

where in the last step the narrow-width approximation 
has been apphed. Using the relation between the pim 
coupling and the partial width 



67r 



as well as the KSFR relation iF^gf, 



indeed 



M2 [2c 



(B4) 



we find 



127r2 



4M? 



e o 



(B5) 



We can also test the sum rule with more realistic in- 
put. As an example, we use the partial wave /i(s) 
based on the finite- matching-point solution ((22|) . with 
C2 = 1.066 X Fstt based on the chiral prediction with 
resonance saturation ([37|) and Sm — (1.2 GeV)^. We 
parameterize the pion vector form factor with a simple 
twice-subtracted Omnes representation 



F,^(s)=exp. 



(r 



,2\V 



s 



ds' Slis') 



(B6) 



4M2 



varying the pion charge radius in the range {r'^)X = 
0.435 . . . 0.450 fm^. Integrating the sum rule (jBip up to 
the matching point Sm saturates the full value F-^^^jT, 
at 87% . . . 90%, depending on the charge radius. The re- 
mainder has to be attributed both to contributions to 
the integral above Sm and to heavier intermediate states 
beyond tttt. Given, however, the degree to which these ef- 
fects (higher energies and heavier states) are suppressed 
already in the unsubtracted dispersion relation, we are 
optimistic that they should be even less important in the 
subtracted form suggested to study the s dependence of 
the transition form factor in Eq. 



In order to obtain a version of Eq. (|13l) that is valid 
off-shell we would have to allow for an additional term 



(s + t + u-3M2)g(s,t,M) 



(CI) 



with an a 'priori unknown function Q{s,t,u). Match- 
ing the twice-subtracted version with one-loop ChPT, we 
find 

C2 + {s + t + u~ Mll)g{s, t, u) ^ F3^{1 + C{s + t + u)). 

(C2) 

This shows that up to higher chiral orders g{s, t,u) = G 
is constant, and that, by comparing coefficients. 



(C3) 



which, unfortunately, merely reproduces Eq. p6l) . 

The same calculation for the once-subtracted disper- 
sion relation, now with a different function Q'{s,t,u), 
gives 



1 



TT Jamb 



d/^^ + g' = F,.c, 



„/2 



(C4) 



and thus again recovers Eq. (p6| . Moreover, the second 
equation corresponds to the sum rule for C from Ref. Q 
if g' is neglected (leading to C = 0.93 GeV"^ Q). Nu- 
merically, this value is rather close to [3] 

^-2M|"96^^-^^°SM|j^^-^^«'^^"^' 

which, however, simply indicates that the off-shell effects 
parameterized by g' are moderate. In particular, the 
HLS expression (|B2)) 



lmh{s)^F3.^M'pS{s-M^p) 



reproduces resonance saturation for C2(/i) 



3e 1 



,Im/i(s') 



3e 



4M2 



1287r2M2 ■ 



(C6) 



(C7) 



Appendix C: Off-shell terms and dispersion relations 



The dispersion relations ((13)) were derived in Ref. [8| 
starting from fixed-i dispersion relations, fixing the sub- 
traction terms by means of crossing symmetry (similarly 
to the derivation of tttt Roy equations [35|). and neglect- 
ing the imaginary parts of partial waves with ^ > 3. 
Although this derivation relies on on-shell kinematics 
s + t + u = 3M^, the same dispersion relations were used 
later on to determine the quark-mass renormalization of 
the chiral anomaly. We will now demonstrate that this 
is not legitimate. 



Besides the issues with the dispersive representation 
for off-shell kinematics, the off-shell ChPT amplitude it- 
self will in general depend on the parameterization cho- 
sen for the pion fields. However, for the present purpose, 
such effects can simply be absorbed into the definition 
of g and g' , whose cancellation in the final expressions 
therefore ensures that Eq. (j36p remains reparameteriza- 
tion invariant. 

All these subtleties regarding the off-shell matching be- 
tween dispersive and chiral amplitudes lead us to the 
conclusion that the quark-mass renormalization cannot 
simply be derived from the sum rules in Ref. (81]. 
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